Search for CP Violation in Charged-S and A Hyperon Decays 
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We have compared the p and p angular distributions in 117 million H~ ^ Att" — * pTT~n~ and 
41 million H"^ — > Avr"*" pn'^Tv^ decays using a subset of the data from the HyperCP experiment 
(E871) at Fermilab. We find no evidence of CP violation, with the direct- CP- violating parameter 
Asa = (ohoa -asaA)/{asaA -I- ohQa) = [0.0±5.1(stat)±4.4(syst)] x 10"". 

PACS numbers: 11.30.Er, 13.30.Eg, 14.20.Jn 



In the standard model, CP asymmetries are expected 
to be ubiquitous in weak interaction processes, albeit of- 
ten vanishingly small. To date, CP asymmetries have 
been seen only in the decays of |l| and Bd mesons 
0. Although the asymmetries observed in these de- 
cays are consistent vifith standard- model predictions, ex- 
otic sources of CP violation have not been ruled out. 
Hence it is vital to search for novel sources of CP vio- 
lation. Hyperon decays offer promising possibilities for 
such searches as they are sensitive to sources of CP vi- 
olation that, for example, neutral kaon decays are not 
0, Q| . The most experimentally accessible CP-violating 
signature in spin-^ hyperon decays is the difference be- 
tween hyperon and anti-hyperon decay distributions in 
their parity- violating two-body weak decays. In such de- 
cays the angular distribution of the daughter baryon is 
dN/dVt = A'o(l + aPp •P(j)/47r, where Pp is the parent 
polarization, pd is the daughter baryon direction, and 
a = 2Re(S'*P)/(|^/ + |^'P), with S and P the ^ = 
(parity-odd) and / = 1 (parity-even) final-state ampli- 
tudes. CP invariance requires that a — —a 0|. 

In HyperCP, the S~ and S+'s were produced at an 
average angle of 0° so that their polarization was zero. 
The angular distribution of p's from unpolarized S~'s in 
S" — * Avr^ — > pTT^TT^ decays is given by 
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since the daughter A is produced in a helicity state with 
polarization as The polar angle 9 is measured in that 



A rest frame, called the lambda helicity frame, in which 
the direction of the A in the rest frame defines the po- 
lar axis. The angular distribution of the p from the corre- 
sponding decay sequence, — > Att^ pK^n^ , should 
be identical if CP is not violated, as both as and aA 
reverse sign. Any difference in the angular distributions 
is evidence of CP violation in either S or A decays, or 
perhaps both. The measured CP-violating observable is 



CtsaA + OiECtA 



(2) 



where = (as -I- as)/ (as — as) and A a = 

(aA + aA)/(aA -oIa)- 

The most recent standard-model calculation for the 
combined asymmetry is -0.5x10~^<Aha<0.5x10~'' 
Note that this prediction uses a theoretical calculation of 
the S- and P-wave Att final-state scattering phase-shift 
differences rather than more recent measurements P]. 
Non-standard-model calculations, such as left-right sym- 
metric models and supersymmetric models |lOl fll| , 
allow for much larger asymmetries. The supersymmetric 
calculation of He et al. J^] generates values of Aa as large 
as 19x10^'*. Bounds from e and e'/e in decays limit 
Aea to be less than 97xl0~^ Q. Experiments have yet 
to probe hyperon CP asymmetries beyond the O(10~^) 
level, with the best hmit being A^a = +0.012±0.014 
|l2l |. In this Letter we present an experimental search 
with significantly improved sensitivity. 

Data were taken at Fermilab using a high-rate spec- 
trometer (Fig.^ ■ The hyperons were produced by an 
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800 GeV/c proton beam incident at 0° on a 2x2 mm^ Cu 
target. Immediately after the target was a 6.096 m long 
curved collimator embedded in a dipole magnet ("hy- 
peron magnet"). Charged particles following the cen- 
tral orbit of the collimator exited upward at 19.51 mrad 
to the incident proton beam direction with a momen- 
tum of 157 GcV/c. Following the collimator was a 13 m 
long evacuated pipe ("vacuum decay region"). The mo- 
menta of charged particles were measured using nine mul- 
tiwire proportional chambers (MWPCs), four in front 
and five behind two dipole magnets ("analyzing mag- 
nets"). At the rear of the spectrometer were two scintil- 
lator hodoscopes used in the trigger: one, the same-sign 
(SS) hodoscope, situated to the beam-left of the charged 
secondary beam, the other, the opposite-sign (OS) ho- 
doscope, situated to beam-right. A hadronic calorimeter 
was used to trigger on the energy of the p or p. 
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FIG. 1: Plan view of the HyperCP spectrometer. 

The (negative) and S+ (positive) data were taken 
alternately, typically with three positive runs followed by 
one negative run in a sequence that usually took about 
12 hours. To switch from one running mode to the other, 
the polarities of the hyperon and analyzing magnets 
were reversed and the targets were interchanged; differ- 
ing target lengths were used to keep the secondary-beam 
rates approximately equal. At a nominal primary proton 
beam rate of 7.5x10^ s^^ the secondary-beam rate was 
13x 10^ s~^, with the average difference between the pos- 
itive and negative rates less than 5%. A simple trigger 
— the "cascade" (CAS) trigger — with large acceptance 
and single-bucket (18.9 ns) time resolution was used to 
select events with the A — > p7r~ topology. It required the 
coincidence of at least one hit in each of the SS and OS 
hodoscopes — the "left-right" (LR) subtrigger — along 
with at least » 40 GeV energy deposited in the hadronic 
calorimeter, an amount well below that of the lowest en- 
ergy p or p. 

A total of 90 billion CAS triggers were recorded in the 

1999-2000 running period. This analysis used data taken 
from the end of the run: a 21-day period in December, 



1999 and a 12-day period in January, 2000. (The inter- 
vening period was devoted to special polarized S runs.) 
The dataset included 19% of all the good S+ events (41.4 
million) and 14% of all the good S~ events (117.3 mil- 
lion) taken in the 1999-2000 running period. The data 
were divided into 18 analysis sets of roughly equal size, 
each containing at least three positive and one negative 
run taken closely spaced in time. 

The data were analyzed by a computer program 
that reconstructed tracks and determined particle mo- 
menta, invariant masses, and decay vertices, assuming 
the S — > Att and A pn hypotheses. Efficiencies of each 
MWPC wire and hodoscope counter were measured on 
a run- by-run basis using tracks from reconstructed 5=*=, 
K^, and fl^ events. These efficiencies were typically 
« 99% and > 99% for the MWPCs and hodoscopes, re- 
spectively. The calorimeter trigger efficiency, as deter- 
mined on a run-by-run basis using good E"^ events from 
the LR trigger, was >99%. Runs with anomalously low 
(<95%) hodoscope, wire chamber, or calorimeter effi- 
ciencies were not used; these were less than 5% of the 
total. The criteria used to select the final event samples 
were: (1) that the pn and pmr invariant masses be, re- 
spectively, within ±5.6 McV/c^ (3.5 tr) and ±3.5 MeV/c^ 
(3.5 cr) of the mean values of the S and A masses (1.3220 
and 1.1158 GeV/c^); (2) that the z coordinate of the S 
and A decay vertices lie within the vacuum decay region 
and that the A decay vertex precede the !B decay ver- 
tex by no more than 0.50 m; (3) that the reconstructed 
S trajectory trace back to within ±2.45 mm (3.3 cr) and 
±3.26 mm (3.4 cr), respectively, in x and y, of the center of 
the target; (4) that the S trajectory trace back to within 
+8.2/ — 8.4 mm and ±6.5 mm, respectively in x and y, 
from the center of the exit of the collimator; and (5) that 
the n^n^n^ invariant mass be greater than 0.5GeV/c^ 
(to remove — > tt^tt^tt^ decays). Cuts on the parti- 
cle momenta and the numbers of SS and OS hodoscope 
hits were also made. Events satisfying these criteria had 
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FIG. 2: The unweighted />7r~7r~ (histogram) and p7r"'"7r"'" (cir- 
cles) invariant masses. 
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a background to signal ratio of (0.43±0.03)% for the ^ 
data and (0.41±0.03)% for the S+ data (Fig. [3). 

The CP asymmetry Aha was extracted by comparing 
the p and p cos 6 distributions in the lambda helicity 
frame. As care was taken to exactly reverse the hyperon 
and analyzing magnetic fields — the fractional differ- 
ence between the magnitudes of positive and negative 
analyzing magnet fields was Ri3xl0~^ — biases due to 
spatial acceptance differences were minimal. The mag- 
netic field magnitudes were updated on a spill-by-spill 
basis using values recorded by Hall probes placed in each 
magnet. Differences in the MWPC wire efficiencies were 
typically on the order of 1x10"'^ in the secondary beam 
region, and much less outside. Hodoscope counter effi- 
ciency differences were typically much less than 1x10^"^. 
These efficiency differences had negligible effects on Aha- 
The calorimeter efficiency difference was « 1x10^"^, and, 
within errors, uniform over the calorimeter face. 

To eliminate differences in the S~ and S"*" momentum 
and position distributions, the and events were 
weighted in the three momentum-dependent parameters 
of the S's at the collimator exit (their effective production 
point): the momentum (pn), the y coordinate (yn), and 
the y slope (sEy). Each parameter was binned in 100 bins 
for a total of 10^ bins. The ps, ys, and ssy bm widths 
were, respectively, 2.25 GeV/c, 0.13 mm, and 0.08x10^'^. 
Bins with fewer than four events of either polarity had 
their weights set to zero. After the weights were com- 
puted the p (or p) cos 9 of each event was weighted ap- 
propriately and the ratio of the weighted p and p cos 9 
distributions was then formed. The expected ratio. 



R = C 



1 + asctA cos 9 
1 -I- (asOA — S) cos 6*' 



(3) 



determined using Eq. was fit to the data to extract 
the asymmetry 6 = a^QiA — asaA — 2asaA-AsA and the 
scale factor C, where a^aA = —0.294 (lj| was used. No 
acceptance or efficiency corrections were made. 

Figure |31 shows typical ratios of weighted and un- 
weighted p to p cos 9 distributions from analysis set 
1. Figure 01 shows S for all 18 analysis sets. Fits 
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FIG. 3; Ratios of p to p cos 9 distributions from analysis set 1, 
both unweighted (filled circles) and weighted (open triangles), 
with fits to the form given in Eq. 0. 



to R were good; the average x^/df was 0.96. The 
weighted average of S for all 18 analysis sets is J = 
(— 1.3±3.0)xl0~'', where the error is statistical, with 
= 24. The corresponding raw asymmetry is 
AHA(raw) = (2.2±5.1)xl0"'^. 

The background-corrected asymmetry was determined 
as follows. The asymmetries in the mass sidebands 1.290 
- 1.310 GeV/c^ and 1.334 - 1.354 GeV/c^ were found, 
using weights from the central region. The weighted av- 
erage of the two sideband asymmetries, scaled by the av- 
erage background fraction of 0.42%, was subtracted from 
the raw asymmetry to give Asa = (0.0±5.1) x lO"**. 

The analysis algorithm and its implementation were 
verified by a simulation, called the collimator hybrid 
Monte Carlo (CHMC), that used momenta and positions 
at the collimator exit from real S~ and S+ events as 
input to computer-generated S decays. Using zero and 
near-zero input asymmetries, the extracted values of Aha 
differed from the input values by (-1.9±1.6) xlO"''. Note 
that the measurement of Asa has no Monte Carlo depen- 
dence. 

Systematic errors were small for several reasons. First, 
common biases were suppressed by taking the ratio of the 
p and p cos 9 distributions. Second, overall efficiency dif- 
ferences do not cause a bias, only spatially dependent 
differences. Finally, since the polar axis changes from 
event to event in the lambda helicity frame, there is only 
a weak correlation between any particular region of the 
apparatus and 9, minimizing biases due to localized dif- 
ferences in detector efficiencies. Tabled lists the system- 
atic errors; added in quadrature they give 4.4x10^**. 

The largest systematic error was due to the uncertainty 
in the calibrations of the Hall probes in the analyzing 
magnets. The magnetic fields were quite stable: varia- 
tions in the Hall probe readings of the sum of the two 
fields were 6.3 G and 5.7 G (rms), respectively, for the 
9024 positive and 2396 negative-polarity spills used in 
this analysis. From calibrations before and after the run- 
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FIG. 4: Asymmetry S (top) and x^/df (bottom) versus Anal- 
ysis Set. 
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TABLE I: Systematic errors. 
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Analyzing magnets field uncertainties 


2.4 


Calorimeter inefficiency uncertainty 




Validation of analysis code 


1.9 


Collimator exit x slope cut 


1.4 


Collimator exit x position cut 


1.2 


MWPC inefficiency uncertainty 


1.0 


Hodoscope inefficiency uncertainty 


0.3 


Particle/antiparticle interaction diflterences 


0.9 


Momentum bin size 


0.4 


Background subtraction uncertainty 


0.3 


Error on as a a 


0.03 



ning periods using more precise NMR probes, the uncer- 
tainty was estimated to be 5.5 G for the sum of the fields, 
corresponding to an uncertainty of 2.4x10"^ in Asa- 

The differences in efficiencies of the calorimeter, ho- 
doscopes, and MWPCs between positive and negative 
running, were not corrected for, as they were negligibly 
small. The effect of calorimeter inefficiency differences 
was determined using a data sample taken with the LR 
trigger. The difference in As.a, with and without the 
calorimeter trigger requirement, was found to be con- 
sistent with zero, with a statistical error of 2.1x10^^. 
Weighting events to correct for the hodoscope ineffi- 
ciencies changed ^ha by only 0.3x10^'*. The effect of 
MWPC inefficiency differences (1.0x10""') was estimated 
using CHMC data by determining the difference in Asa 
using real and 100% efficiencies. 

The effect on A-^a of tighter cuts on the (unweighted) 
S X slope and position at the collimator exit was stud- 
ied and resulted in respective uncertainties of 1.4x10"'' 
and 1.2x10"''. The effect of the bin sizes used in ex- 
tracting the event weights was investigated by increasing 
and decreasing the S momentum bin sizes by 25%, A^a 
being most sensitive to momentum. Another possible 
source of bias was a momentum-dependent differential 
loss of events due to interactions of the S~ and S+ de- 
cay products with material in the spectrometer. Monte 
Carlo studies, using the interaction cross sections given 
in Ref. 0|, showed this bias was negligible. 

The result was stable with respect to time, S mo- 
mentum, and secondary-beam intensity. Differences be- 
tween the S~ and S+ production angles could in prin- 



ciple cause a bias due to production polarization differ- 
ences. Average production angle differences were only 
w 0.02 mrad. Assuming a linear dependence of the polar- 
ization on transverse momentum |l5l | , Monte Carlo stud- 
ies indicated a negligible effect on the p and p cos 9 slopes. 
No dependence of A^a on production angle or incident 
proton beam position was evident. 

To conclude, we have measured Asa to be 
[0.0±5.1(stat)±4.4(syst)] xlO"''. This result is consistent 
with standard-model predictions and is a factor of 20 im- 
provement over the best previous result [l^. 
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